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Exact solutions of the Schrodinger equation for the Coulomb potential are used in the scope 
of both stationary and time-dependent scattering theories in order to find the parameters which 
define regularization of the Rutherford cross-section when the scattering angle tends to zero but 
the distance r from the center remains fixed. Angular distribution of the particles scattered in the 
Coulomb field is investigated on the rather large but finite distance r from the center. It is shown 
that the standard asymptotic representation of the wave functions is not available in the case when 
small scattering angles are considered. Unitary property of the scattering matrix is analyzed and 
the "optical" theorem for this case is discussed. The total and transport cross-sections for scattering 
of the particle by the Coulomb center proved to be finite values and are calculated in the analytical 
form. It is shown that the considered effects can be essential for the observed characteristics of the 
transport processes in semiconductors which are defined by the electron and hole scattering in the 
fields of the charged impurity centers. 
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I. INTRODUCTION 

Scattering of non-relativistic charged particles by the Coulomb center is one of the canonical problems both in 
classical and quantum mechanics which is known as the Rutherford problem. It is the standard point of view that 
the differential cross-section da{6) of the particle scattering to the solid angle dVl has the same form in the both cases 
(for example, Refs. 0,0) 

da{e) ^ a{6)dn = {^f-^. (1) 
2mv^ sm^/2 

Here m and v are the particle mass and velocity correspondingly, parameter a defines the amplitude of the Coulomb 
potential U{r) — ajr. 

So, the main measured characteristic of the scattering process in the Coulomb field has the non-integrable singularity 
in the limit 6* ^ (in quantum theory the singularity exists also in the scattering amplitude). Fortunately, this 
singularity doesn't lead to any problem when describing of the most real experiments because particles are scattered 
by the systems with zero total charge. In this case the singularities conditioned by the scattering centers of opposite 
signs are compensated and the cross-section proves to be regular in the entire angular range. Nevertheless, there 
are some physical systems where one should consider the problem of regularization when calculating such integral 
scattering characteristics as the total atot and transport Otr cross-sections 



otot = j da{9), atr = cos 9)da{9). (2) 

As for example, we can mention calculation of the characteristics of kinetic processes in plasma and impurity semi- 
conductors or collisions of the charged particles in beams. In such cases one should introduce some phenomenological 
parameter Omin for cutting off the cross-section with angles 9 < 9 mm- This parameter can be defined by various 
physical reasons. Particularly, in the framework of the classical mechanics the small angle scattering is defined by the 
particles with large impact parameter 3] connected with a long range character of the Coulomb potential. Therefore, 
the small angle cone can be excluded from the consideration because of finite transversal width a of the incident beam 
with 9min a/r 0]. 

Another approaches are used when the mobility of the charge carriers is calculated in the impurity semiconductors. 
The models of Brooks-Herring and Conwell-Weisskopf are mostly used for this problem at present. These models 
correspond to different ways for estimation of the parameter 9min connected with screening of the Coulomb potential. 
However, such estimations have only qualitative character and some additional phenomenological parameter should be 
introduced for more precise description of the mobility as it was shown recently in the paper . Accurate calculation 
of the integral values characterizing the charge carrier scattering by impurities is actual because of high accuracy of 
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measurement of these values in real semiconductors (for example, |^). Solution of this problem is of great interest 
also for analysis of the electron transport in nanostructures such as quantum wires 0], superlattices and films |lClj |. 
nanotubes [U^ . 

Regularization problem for the Coulomb cross-section is essentially more principal in the framework of the quantum 
theory. The matter is that the exact wave function for the states of the continuous spectrum is well known Ij and 
it has no any singularity even in the case of the plane incident wave which corresponds to the beam with the infinite 
transversal width. It should mean that the singularity of the scattering amplitude is not intrinsic feature of the 
Coulomb system in the scope of quantum mechanical description. Possibly, it could be conditioned by not completely 
adequate interpretation of the asymptotic behavior of the wave function in this case. One can expect that some 
characteristic, or "kinematical" , regularization parameter 9q should exist which doesn't connect with the initial state 
of the system unlike the value 9min- In general case the regularized cross-section should depend on both parameters. 

It is essentially to emphasize that some specific characteristics of the Coulomb scattering problem have been widely 
discussed in monographs and textbooks. As for example, it was shown in the book jj] that the connection between the 
impact parameter and scattering angle becomes indefinite in the case of 6* = 0, therefore the scattering cross-section 
for zero angle can't be calculated in classical dynamics. It is also well known that long-range character of the Coulomb 
potential leads to the logarithmic distortion of phase in the asymptotic form of the wave function (for example 1]). 
However, the problem of the cross-section regularization has not been considered in these discussions. 

This question was analyzed for the first time in our paper |l2j . It was shown that the standard asymptotic 
representation of the wave function was not actually formed in the range of small angles when considering the scattering 
processes by the long-range potentials {U{r) ~ l/r'*;s < 3). In the result the canonical definition of the scattering 
amplitude proved to be unavailable. Born approximation over the potential U (r) and the non-stationary collision 
theory j^^l were used in our wor k in order to calculate the scattering cross-section without any singularities. We 
can also mention several papers ( and references therein) where it was shown that the interference between incident 
and scattered waves changed the asymptotic form of the wave function and could be essential in real experimental 
conditions even in the case of some short-range potentials. 

In the present paper we consider the non-asymptotic analysis of the observed characteristics for the non-relativistic 
Coulomb scattering problem out of the framework of the perturbation theory. We use the exact solutions of the 
Schrodinger equation in order to answer the following questions: l)which "intrinsic" kinematical parameter defines 
regularization of the Rutherford cross-section in the framework of the stationary scattering theory; 2) how does 
this regularization depend on "external" parameters such as the transversal width of the incidence wave packet or 
effective cutting off of the potential; 3) which way can one calculate non-asymptotic values for the integral scattering 
characteristics crtot,crtr; 4) what is the analog of the "optical" theorem 47r/m/(0) = katot Q, the case of 
the Coulomb potential? It seems to us that the answers for these questions have the important methodical value 
for understanding the scattering processes in the field of long-range potentials but have not been discussed earlier. 
Besides, these results can be also essential for some applications such as the above-mentioned transport processes in 
the semiconductors with the charged impurities. 

The paper is organized as follows. In Sec. 2 the differential scattering cross-section is defined without asymptotic 
representation of the wave functions and the kinematical regularization parameter is found for the Rutherford problem. 
The most important integral characteristics of the scattering problem are calculated in Sec. 3. In Sec. 4 the scattering 
operator and the conservation of the total flux are analyzed. The time-dependent consideration of the collision process 
is discussed in Sec. 5 and influence of the incident beam parameters and screening of the potential to the observed 
scattering characteristics is estimated. The scattering characteristics of the carriers in non-degenerated semiconductors 
with the charged impurities are calculated in Sec. 6 and the results are compared with the experimental values of the 
carrier mobility in real systems. 



II. NON-ASYMPTOTIC CALCULATION OF THE DIFFERENTIAL CROSS-SECTION FOR THE 

COULOMB SCATTERING 

Let us remind the standard deflnitions of the scattering theory in the stationary quantum mechanics. It is well 
known I*], that in this case the wave functions of the continuous spectrum iljj:{r) should be found as the solutions of 
the Schrodinger equation 

[-|l(A + fc2) + t/(f)]^_(F)=0, (3) 
with the following asymptotic boundary conditions (Fig.l shows all necessary notations). 



FIG. 1: Sketch of the scattering process in the stationary case. 



V'fe(^)~e"=^ fcf^-oo; (4) 



V'fe(^1~e^"'' + /(e) — , r»i?. (5) 

Here k is the wave vector; the value R defines characteristic radius of the potential action with the center point 
r = (i? ^ oo in the case of the Coulomb field); the wave function is supposed to be normalized to one particle, so 
that the flux density in the incident state is: 

J=^[^fc(^*Wfe(?=0-^fe(^Vi/.j(f)*]~^^ = Jo, fcf^-oo; (6) 
The flux density in the asymptotic state ((SJ is divided on two components: 

/-i/7+J.c-, r>i?. (7) 
k r 

One of them ji (longitudinal component) corresponds to the particles passed through the field without interaction 
and the second one jsc (radial component) describes the scattered particles. It leads to the standard definition of the 
cross-section: 

da(e) = —r^dn; 
Jo 

^{0) = \f{0)?. (8) 

It should be noted that the longitudinal flux is also changed ii < jo, and its decrease is deflned by the total 
scattering cross-section in accordance with the " optical" theorem [l| . 

Evidently, that the deflnition ^ is based essentially on the asymptotic regime (jSJ for the wave function in the 
observation point r. Accordingly to the terminology used in radio-physics and optics (for example, |l5j|'). it means 
that the particle should go out of the "near" zone, where the action of the potential is still essential, and pass to the 
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"far", or "wave", zone. The boundary between these zones is defined by the condition that the interference between 
incident and scattered waves becomes neghgible that is the difference between their phases satisfies the inequahty 



kr - kf^ 2krsm'^{9/2) > 1; 

0>eo = J^. (9) 
V kr 

We suppose further that for aU real coUisions the condition fcr 1 is fulfilled. 

It is clear that the boundary of the " wave" zone depends both on the distance r from the center and the scattering 
angle 9 (Fig.l). It means that in general case there is a part of the particle fiux which can not be described by the 
asymptotic wave function © even for rather large distance r. Certainly, that this property does not depend on the 
radius of the potential auction. However, the question is: what is the contribution of these particles to the integral 
scattering process? When the distance from the center r is fixed, the number of particles scattered to the "near" zone 
6 < Oq can be estimated as 



iV,,,c.jo-(O)0„2-^, (10) 

The cross-section a{0) is restricted for the potentials with the finite action radius R, therefore the Ndif decreases 
quickly at the large distance. It means that the contribution of these particles to the observed scattering characteristics 
is negligible for the most real experiments. The detailed analysis of the "near" and "wave" zone formation for the 
scattering problem with the short-range potential has been recently considered in the paper |l4j . 

The picture changes fundamentally in the case of the long-range potential (R — > oo) . The value NiHf can even 
increase with the distance and its contribution to the formation of the scattering flux can be essential. Particularly, 
the analogous estimation in the case of the Coulomb field leads 



^-/»^o(^f^-^.. (11) 

It means that the asymptotic boundary condition ((Sj is not available in the entire range of the scattering angles 
and the nonasymptotic expression for the wave function should be used in the case of small angles. It is important 
to stress that this circumstance does not connect with the width of incident beam and defines by the characteristic 
feature of the potential itself. 

So, the considered regularization problem for the Rutherford cross-section in the scope of the stationary scattering 
theory is reduced to the analysis of the space fiux distribution on the basis of the well known exact solution of the 
equation with the potential U{r) — Ze^/r but without handling to the asymptotic representation of the wave 
function. 

We will use the following form of the normalized wave function Q 



^P^{f) = Ne'''^F[±i£„ l,i{kr - kf)]; N = e^^«r(l =F 

^ a „ o hk , . 

a = Ze^; v = —, 12 
hv m 

where F(a,b,t) is the confiucnt hypergeometric function; T[t) is the Gamma - function; the upper sign in the 
formulas corresponds to the attraction field and the lower one corresponds to the repulsion potential. 

Let us show that the fiux density in the formula ^ calculated with the exact wave function can also be divided 
by two components in accordance with the formula ((TJ as it was in the asymptotic regime . For this purpose one can 
use the following representation of the function F as the superposition of two confiuent hypergeometric functions of 
the 3-rd genus J7i_2(a,6, t) 0| 
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Fi±it I, iz) = -I_C/i(±z^, 1, ^z) + —1-—U2{±tt 1, iz); 

r(±«^) r(i=FO 



Jo i-z 



du 



Jo 

z = {kr — kr) = kr{l — cos 9). (13) 

Let us also mention the connection between these functions and the confluent hypergeomctric functions of the 2-nd 
genus U{a, b, t) ^ 

Ui{a, b, t) = U{b - a, b, -i)e*e±"(°-''), 

1 

U{a,b,t)^—— e~*''u''-\l + u)''-''-^du;. (14) 



r(«) Jo 

When the Rutherford cross-section is calculated by means of the standard definition this representation permits 
one to find the asymptotic form of the wave function in the hmit 2; ^ 1 This case corresponds to the "wave" 
zone when the function Ui transforms to the spherical wave and the function U2 tends to the plane wave. However, 
both these functions are well defined also in the "near" zone (z < 1) when they can be calculated by means of the 
following serieses p^: 

Ui{±it 1, iz) = ^T{±iO{F{±iC, 1, iz) + 
[[21n(iz) T «7rcoth«) - in + 2^(±i£,)]F(±i£,, 1, iz) + 

m=l ^ ^'^ 

U2{±it 1, iz) = ^r(l T ^0{F{±iC, 1, ^z) - 

[[21n(iz) =F i7rcoth«) - in + 2i:(±i£,)]F(±iS,, 1, iz) + 

2ni 

^ £ nlTeKmO^ t^^"^ ^ " "^^^'^^ + ^^(^^ " + m^z)^}, (15) 

where the logarithmic derivative of the Gamma-function. 

When the representation H13|l is used in the formula © one should take into account only the derivatives from the 
exponents because the conditions fcr 3> 1; z ^ 1 supposed to be fulfilled. As for example, 

-*(-V)[e''^''-Gi(zz)] ^ e^^'-ifcGi + -G'Mz)] ~ e*'='-fcGi[l + 0(-^-)]. 
r r Kr 

This representation permits one to find the scattering fiux jsc directed to the observation point along the vector r 
without use of the asymptotic form ([SJl of the wave function. In the result the scattering cross-section can be defined 
in the entire range of the angles 9 in the following form: 
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ai{d)de sin 9r^ / d^—de = 
Jo jo 

2i sinh(7rOe=F'^« |Gi (ii^ , iz) jV^ sin OdO. (16) 

One can see that the differential cross-section ai{9) is finite for any angle in spite the function Gi(±i^,z2;) has the 
logarithmic singularity at zero angle as it follows from the equation But this value depends on the distance 

between the center and observation point by non-trivial way because of long-range action of the potential field to the 
particle. The result of this action at small angles ("near" zone) does not reduce to varying the phase of the scattering 
amplitude as it is takes place for the asymptotic range of angles ("wave" zone)0|. 

If one considers behavior of the function G'i(±i^,zz) in dependence on the scattering angle the " kinematical" 
parameter 9q for regularization of the Rutherford cross-section can be introduced by the natural way. Actually, the 
asymptotic range of angles corresponding to the "wave" zone is defined by the condition 



1 2 

z = fcr — fcr ~ — fcr^ ^ 1; 
V kr 



—; z~x. (17) 

Here the dimensionless value x is introduced as the convenient variable for the angles compared with the width 
of the "near" zone. Certainly, in the range of a; 3> 1 the standard asymptotic representation of the integral in the 
definition of the function Gi(±i^, iz) leads to the result corresponding to the formula with new variable 



ai(g)^27r(|)^ ^^y^' ^8^(^)^l; (18) 

It is well known that the interference between the scattering flux and the flux directed along the initial velocity 
of the particle does not take into account in scope of any quantum scattering theory based on the solutions of the 
stationary Schrodinger equation [H. So, in order to use the formula (|16|l in the range x < 1 corresponding to the 
"near" zone, one should compare it with the angle width of the zone where the above mentioned interference is 
still essential. It is clear that the angle width of such "interference" zone does not depend on the dynamics of the 
interaction between the particle and field. It is defined only by the transversal width a of the incident particle wave 
packet (Fig.l). One can estimate the angle width 9int of the "interference" zone as follows 



0^nt ^ -. (19) 

r 

It means that one can consider the scattering flux in the near zone and at the same time neglect by its interference 
with the incident beam if the following conditions arc fulfilled 



dint < 9 < 9o — 




These inequalities are satisfied in the case of rather large r as it usually supposed in the scattering theory. More 
accurate analysis of this factor will be considered below (Sec. 3) in the framework of the time-dependent theory of 
collisions But one can estimate just now the contribution of the "interference" zone to the integral scattering 
characteristics which are finite values in our consideration unlike the asymptotic analysis. As for example, the ratio 
of the particle flux scattered at the "interference" and "near" zones can be estimated as 



S=^^ ai(0)d0/ / " ai{9)d9 =^ ^ « 1. (21) 

Jdif Ja Jo 2r 

It remains small under standard conditions of the collision theory [T^ and one can analyze distribution of the flux 
density in the "near" zone neglecting its interference with the incident flux. It permits one to find the leading terms 
of the differential scattering cross-section at small angles using the series (|15|) 



7 



(71 (61) ~ 8\/2^e=F'^« sinh(7rOa;(lna;)2^. (22) 

V A: 

Fig. 2 compares the accurate and asymptotic scattered fluxes for various values of the variable x and parameter ^. It 
is interesting to pay one's attention to the essentially different behavior of the nonasymptotic flux in "near" zone for 
scattering by the attractive and repulsive centers in contrast to the Rutherford cross-section ^ which is independent 
of the potential sign for any value One can see that the regularized differential cross-section (|16() in the "near" 
zone is essentially non-invariant relatively to the sign of the charge if the parameter ^ > 1. It should be noted that 
the effect of slightly different interaction of the charge carriers with the impurities of different signs is well known 
in the semiconductor physics. It is usually considered there by means of the Friedel sum rule using the partial 
expansion of scattering amplitude in the series of orbital momenta. 

As it follows from Eq. 1)22(1 , the scattering flux in the case of attractive potential varies rather slowly with increase 
of the parameter ^, but it grows exponentially in the case of repulsion. Certainly, such behavior of the cross-section 
takes place only in the narrow angle domain H2UII and compensates the exponential decrease of the flux just along the 
line 6 = {) which is well known for the repulsive potential [l| 
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FIG. 2: Ratio of the scattered flux J„ to the density of the incident flux Jo. Solid line - the case of attraction Jatt, dashed 
line - the case of repulsion Jrep, dotted line - the Rutherford flux J ruth- 



III. INTEGRAL CHARACTERISTICS FOR THE COULOMB SCATTERING PROBLEM 

Let us now calculate the integral scattering characteristics for the considered problem. In accordance with Eq. H16|l 
the nonasymptotic expression for the total cross-section is deflned by the following integral 

7r 

2^ sinh«) I Gi(±iC,zz) I V sin 6ld6'. (23) 

One can use in this integral new variable z 
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utot = j C sinh(7r^)|G'i(±i^,iz)p— dz, 



(24) 



and represent it as the sum of two integrals 



CTtot = -^Csinh(7rO{ J \Gi{±i^,iz)\'^dz - 

\Gi{±i^,tz)\^dz}. 



(25) 



2kr 



The asymptotic representation for the function Gi can be used in the whole interval of integration in the 
second integral 



r(iT*0 



and it leads to the following simple result 



2r f°° Tvf'^ 

h = ^esinh(7rO / \Gi{±i£,M?dz ~ 

/C J2kr 



(26) 



Therefore it has the order of {kr)~^ in comparison with the first integral and its contribution to the total cross- 
section can be omitted. It permits one to find how the value atot depends on the most essential parameters of the 
problem 



/±(0=e^"M \U,{l±i^,l,tz)\^dz. 
Jo 



(27) 



Here we use again the canonical form for the confluent hypergeometric function of the 2-nd genus pla|. The universal 
functions I± (f ) depend only on the variable ^. They are defined by the converged integrals and can be easily calculated 
numerically. Fig. 3 shows the results of these calculations. 
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FIG. 3: Universal functions /±(C) which define the total scattering cross-section as the functions of the parameter ^ for the 
cases of repulsion and attraction. 

Now let us consider another integral characteristic of the scattering process, namely, the transport cross-section 
which is very important value for a lot of applications. It is defined by the formula 



atr= f 2^smh{Tr£,)\Gi{±i^,iz)\^r^sm9{l~cos 
Jo 



(28) 
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If one uses the variable z in this integral and comes back to the hypergeometric function of the 2-nd genus, Ea. H28|) 
transforms as follows 



r2kr 



4^(0 = / \U(l±i^,l,tz)\^zdz. (29) 
Jo 

The integrand function for the transport cross-section is essentially suppressed in the range of small angles in com- 
parison with the total cross-section. Therefore Ut,- is defined by only the logarithm of the distance to the observation 
point unlike to atot proportional to this distance. Besides, this function decreases rather slowly for the large z and we 
can't use the trick analogous to Eq. (|25|l for atr- Nevertheless a series of transformations of the integral I±{£,) permits 
one to find the analytical dependence on the coordinate r with an accuracy of the order (fcr)^^. Let us separate the 
integral on two parts by the following way 

4(0=/ \U{l±i£„l,iz)\^zdz+ \U{lTi^A,iz)\^zdz. 
Jo J I 

If one uses the asymptotic formulas for the hypergeometric functions |0 

|C/(l±z01,zz)P~^+O[(A:r)-3], 
z'^ 

the second term in this integrals is transformed identically 



\\U(\±ii,\,iz)Y / dz 



The second integral here is calculated analytically but now the integrand expression in the first one decreases rather 
quickly and the estimation analogous to can be used 



/ [|C/(l±l01,iz)|2 —]zdz~ / [|C/(l±i01,Zz)|2 —]zdz~ 

Jl Z J2kr Z 

[\U{1 ± i^, l,iz)\^ —]zdz + OUkry^]. (30) 

z^ 

In the result the transport cross-section is defined by well converged integrals 

^ -{ / \V{\±ii,\,iz)fzdz + 
Jq 



CTtr 



fc2 



poo ±7ri; 

J [\U{1 ± itl, iz)f - ^]zdz + e'^''^\n{2kr)}. (31) 
Fig. 4 shows the results of numerical calculation of the universal functions I±{^) — [k'^crtr ~ Stt^^ ln(2fcr)]. 

IV. SCATTERING OPERATOR AND CONSERVATION OF THE FLUX IN SCOPE OF THE 

STATIONARY THEORY 

As it follows from the results of the preceding section the integral scattering characteristics calculated on the basis 
of nonasymptotic consideration increase together with distance r to the observation point. It seems for the first sight 
that it can contradict to the conservation of the total flux of the particles when r becomes rather large. However, let 
us show that this dependence expresses only the fact that the potential influences on the scattering process at any 
distance from the center but the scattering flux remains essentially less than integral incidence flux at any r. For 
qualitative analysis one should take into account that in scope of the stationary scattering theory the quantum state 
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FIG. 4: Universal functions I±{^) for the transport cross-section in the cases of attraction and repulsion. Solid line - the case 
of attraction, dashed hue - the case of repulsion 



of the incident particle is described by the plane wave. Then the total incidence flux Jq through the sphere with 
radius r corresponding to the observation point can be estimated as follows 

Jo - joTrr^- 

Then the ratio of the scattering and incident integral fluxes is 



■tot o 

j^^at^^ /± « 1. (32) 

It is also important to consider this problem more precisely. It is known T| that the condition of the total flux 
conservation leads to the "optical" theorem in the quantum theory of scattering by short-range potential when the 
amplitude of scattering to the zero angle is the finite value. We use the same approach 1] in order to find the 
consequence of this condition in the case of nonasymptotic analysis of the Coulomb scattering. 

Let us represent general solution of the Schrodinger equation in the case of the elastic scattering as the linear 
combination of the functions (|12ll with arbitrary coefficients $(n) which define the amplitudes of probability to find 
the state with the wave vector fc = fcn in the initial packet: 

^{r)= J $(n)V'fe(r)da3 = iV J ^{n)e'''''^" F[±i(,l,ikr{l ~ nn')]dnn = 

pikr ^ikrnn' 

T 

ft = — ; z — fcr(l — nn'), (33) 
r 

dflf} is the element of the solid angle in the direction of the vector n. 

In accordance with the physical interpretation of the contributions defined by the functions Gi , G2 to the total wave 
function p2|l . the first term in Eg . H33|l describes that part of the integral scattering operator 1] which corresponds 
to the formation of the scattering wave. The term, proportional to the function 6*2, describes the deformation of the 
wave packet conditioned by change of the plane wave in the Coulomb field. One can estimate the second term by 
the same method that was used for proof of the " optical" theorem in the case of short-range potential Q . If the 
condition fcr 3> 1 is fulfilled, the main contributions to this integral are defined by the small intervals near the points 
of the stationary phases when integrating over n . These points correspond to the vectors ni = —n' and n2 = n' . Near 
the first point the variable z ~ 2kr is very large. Therefore one can use the asymptotic expression for the function 
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G'2(±i^,iz) and the integrand has no singularities in this case. In the result the contribution to the integral from the 
domain close to this point defines the converged spherical wave with the standard logarithmic distortion of its phase 

Q 

^ — ikr^i$^ In 2kr 
- 2TTi ; $(-n')- 

kr 

When estimating the contribution to the integral from the second point of stationary phase corresponding to 
scattering at small angles one should take into account that the function G2{ii£,,iz) has the logarithmic singularity 
in the point z = 0. Nevertheless, the rather smooth weight function <f>(n) can be removed from the integral in the 
point n — fi' . It leads to the following estimation: 

ikr poo noo 

^2^$(n)eT^«— — — [/ {z)^'^e-''G2{±it^z)dz- / {z)^'^e-'^G2i±i^,iz)dz]. 

The second integral in this expression can be omitted in the limit fcr 3> 1 and the initial wave function is represented 
in the form: 

^{f)~2m <^{-n') - 2TTi-—[A<^{n') + / f{n,n')<^{n)dn]; 

kr kr J 

A = / {z)^'^e-''G2{±ii, iz)dz- 



/("'^') = -^^^eTt«e±^«i"'='-(l~nn')±'«Gi[±z^,ifcr(l-^^')]- (34) 
2ttT[±iQ 

It is more convenient to rewrite this expression in terms of the hypergeometric function of the 2-nd genus 



/•oo 

v4 = ieTf« / e~'''U{±i^,l,iz)dz] 
Jo 



f{n, n') = -!^J^iil^eT^«C/[l T ^t 1, -*A:r(l - nn')]. (35) 

Now the function is represented as the superposition of the ingoing and outgoing spherical waves and it permits one 
to introduce the scattering matrix [? ] as the following integral operator: 

Sin,n')~ASnM'+fin,n'). (36) 

Here 5^,™' is the unit operator which corresponds to the wave passed without scattering and the parameter A 
defines the change of its amplitude (in the case of the short-range potential A = 1 0]). The integral over angles from 
the operator J f*{n,n')f{n',n)dn' coincides exactly with the expression for the total cross section H27fl . Long-range 
character of the potential is appeared in the fact that the scattering matrix elements depend on the coordinate r. 
However, it is very important to introduce such operator because just it defines the kernel of the collision integral 
in the kinetic equations for description of various transport processes But if one uses such operator in the 

collision integral for one-particle distribution function the additional averaging over the coordinate should be fulfilled. 
Dependence of the function f(n',n) on the coordinate is rather smooth , therefore the value r in this function can 
be substituted as an average distance between the scattering centers if the correlation between these centers can 
be neglected (see below §6). Analogous substitution was used in some well-known models for regularization of the 
transport cross section of scattering by the charged impurities in semiconductors , . 

Unitary property of the matrix S{fi,n') leads to the "optical" theorem in the case of short-range potentials 
But if one uses this condition in the case of Coulomb potential there is the problem that the operator f{n, n') has 
the logarithmic singularity in the limit of coinciding arguments and one should define the way for calculating integral 
from the product of singular functions /(n, n') and 5s,n' in the operator SS~^ . Actually it means that the asymptotic 
estimation of the integral in Ea. l|34|) is unavailable for the operator which is quadratic over the scattering matrix. 
Therefore let us analyze separately the conservation of flux considering the following integral 
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I = J dS{Vjir)) = y ((fj(f))rdat; 

= - *(r)V**(f)}, (37) 

with the total wave function H33|l . 

When the superposition H33|) is used in formula H37|l one can take into account the completeness of the coefficients 
$(n). Then integration over all directions in this integral is equivalent to the integral from the flux j^t calculated by 
means of the general formula ^'^7^ but with the stationary wave functions i'j:{r) defined by Eq. I|12|l ( let us consider 
the attractive potential for definiteness ) 



J.t(r) = ^^—Ak\F[i^, l,^{kr~kr)]\' - ^{k - k-)%FF*)}- 
m smh 7r4 r 

VF[i^, 1, i{kr ~ kj^] = - k-)F[i^ + 1, 2, i{kr - k?^] = ^{k ~ k-)Fi. (38) 

r r 

Certainly, the value / is equal to zero identically because of the flux conservation for the stationary scattering 
problem. The "optical" theorem is followed from this condition if the asymptotic form Q for the wave function can 
be used 0. But in the considered problem this condition means that the flux directed along the vector k (it defines 
change of the intensity of the incident wave) , and the scattering flux along the vector f are connected as follows 



j dnfi{kn)\F[iS,,l,i{kr -kT^W^ = i j dnn{kf-k)^{FF^). (39) 

As it was shown above the integrals over the angle for the Coulomb scattering problem include essential contribution 
defined by "near" zone. Therefore the both parts of Ea. H39(l depend on the coordinate r and the standard asymptotic 
expressions for "optical" theorem is inapplicable because the total cross section and the scattering amplitude at zero 
angle are tending to infinity in this case. But if one shows that the leading terms of Ea. (|39(l are equal in the limit of 
large r {kr ^ 1) it can be considered as the analog of the "optical" theorem for the Coulomb potential. 

In order to prove it let us use new variable for the integrals in Eq. 



dz 

kr — kr; sin 9d9 — - — , 
kr 



and transform them as follows 



/ |F(*e,l,iz)pdz- / f {|F|2+e3[F(-*^,l,-*z)F(*e+l,2,zz)]}dz. (40) 
Jo Jo kr 

One can estimate the integrals from the confluent hypergeometric functions in the range fcr 3> 1 by means of the 
following approach. Integral in the left side of Eq. (|40|l can be transformed identically 



p2kr /"Oo poo 

= / \F{i^,l,iz)fdz = \im[ iFfe-^'^dz- iFfe'^^'dz]. (41) 

JO "5^0 Jq J2kr 

The parameter (5 — > is introduced for the regularization of both integrals at upper limit. The asymptotic form of the 
function F can be used in the second term and the first term can be expressed through the hypergeometric function 
F{a, P, 7, z) by means of the formula (see, for example Ij) 

/•oc 

J(A)= / e-^'z''-^F{a,-f,kz)F{a',-f,k'z)dz = 







r(7)A"+" -^(A - fc)-"(A - k')-"- F[a, a', 7, - k)iX - k') ^' ^^^^ 
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When the integrals from the functions with different second arguments are calculated, the following recursion 
relation can be used 

F{a + 1, 7 + 1, z) = + 1, 7, z) - F{a, 7, z)]. 

Let us write also the leading terms of the asymptotic expansions for the functions F Fi which are used for the 
integrals in the limits {2kr, 00) 

FUf 1 ^z) - (1 - ^ + i!ii±^) - !^]- 

F, . + 1, 2, . ^[^.(^(1 - - ^.(^(1 - ^^)]; 

In the result the leading term in the left side of Eq. H4U|) is the following 

fK^"^' ' 'r(-.o 'll + ^lll^l- <«) 

This value defines variation of the flux directed along the incident wave vector and it grows linearly together with 
the distance from the scattering center analogously to the total cross section. As it was mentioned above (Eq. (j2Hl)), 
this growth is not connected with increase of the particle flux but describes distorted part of the wave front which is 
extended together with r because of long-range character of the potential. 

Calculation of the integral 

r2kr 



J2 = — \F{iix%z)\^zdz 

by means of the analogous technique leads to the following result 

TT^ kr 1 



where '^{x) is the logarithmic derivative of F - function p^ . 
The last integral in Eq. (gOJ 



J3 = -L / z3[F(-zC, 1, -iz)F{i^ + 1, 2, iz)]dz = 
«^ Jo 

5R[i^(-i^, 1, ~iz)F{i^ + 1, 1, (45) 



^ ^2fcr 

kr 



transforms as follows 



J3 ^ -f e--«^^{e^(-5R^(l + ^0 + ln2fcr) - 
kr 7r4 

Substitution of Eqs. (|43|l - (|46|l to Eq. H40I) shows that it is satisfied with the considered accuracy. Besides, one 
can see that the left side of Eq. H4U|I coincides with the total cross section (|27|l in the limit fcr ^ 1, and the right side 
of Eq. (|40|) transforms to the imaginary part of the scattering operator H4U|) with n = ft' . So, we can consider this 
calculation as the proof of the "optical" theorem for the Coulomb scattering problem. 
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V. MOVEMENT OF THE WAVE PACKET IN THE COULOMB FIELD 



As it follows from the results of the preceding sections, regularization of the Rutherford cross section is defined by 
the characteristic angle 




which corresponds to the boundary of "near" zone and is considered as the kinematic parameter (KP) of the system. 
However, in real scattering experiments the incident particle is actually represented by the localized wave packet [l3j| . 
Besides, the Coulomb potential is screened at some distance Rg, depending on the properties of the medium where 
the collision is happened. Therefore in general case the problem is characterized by some additional parameters that 
can be considered as the external parameters (EP). So, it is essential to estimate the conditions when the KP is more 
important for the cross section regularization that the EP. We will take into account two the most essential EP: the 
screening angle 9s the incident angle parameter 9i„t , depending on the wave packet transversal width a and defining 
the zone of interference between the incident and scattered waves (see also Sec. 2). The simple estimation of these 
parameters leads to 



1 

kRs 



a 
r 



(48) 



Evidently, the kinematic regularization is the most essential if the angle width of the near zone is larger in comparison 
with the characteristic angle intervals connected with EP, that is the following conditions are fulfilled 



9o>es, ^>l; 9o>e,^t, — <1. (49) 
r r 

The first inequality depends on the mechanism of screening and should be analyzed for every concrete system as it 
will be considered below (Sec. 6) for the scattering by impurities in semiconductors. In order to take into account the 
finite size of the wave packet in the second inequality in (|49|l one should use the time-dependent theory of collisions 
HH; that we will consider in this section. 

Let us suppose that the initial state of the particle in the moment i = is defined by the wave packet in the 
following form 

*-(r,0) = J dq^[q- k)e"^<-''-''°^ = e'^(""~"^°)G(|f - fo|), 

G{p) = / dp^{p)e'PP, (50) 



where fo is the coordinate corresponding to the initial position of the wave packet ; $(p) are the amplitudes of 
probabilities of the wave vector distribution near the center k in the initial state; G{p) is the function which describes 
the form of the localized wave packet in the coordinate space p^ . 

In order to describe evolution of the wave packet (|5U|I it should be expanded in the solutions of the stationary 
Schrodinger equation with the Coulomb potential 12] (let us consider the attractive potential for the definiteness ) 



V'£(r) = Ne'^''F[i^k, 1, «(fcr - kr)]. 

In the standard experimental setting (Fig. 5) the initial position of the wave packet corresponds to the condition 
zq — > —oo. In this case the stationary wave function 'ipj:{r) coincides with the plane wave jlSj and the expansion of 
^i^{r, 0) in the functions ip^{r) includes the same coefficients as in the formula (|5Uf) with an accuracy to the terms of 
the order jzol^^ conditioned by the logarithmic distortion of the wave front in the Coulomb field In the result the 
wave function describing the wave packet state in an arbitrary moment of time has the following form 



^j:{r,t) = / dq<^iq- ^)e'9~e-'«"«e"«'/2r(i _ i^g)F[i^g, l,i{qr - qr}]e-'^^'*/^' 



(51) 
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As it was investigated in detail in the monography [I^ the wave packet spread (diffraction) can be neglected 
during time of the interaction in real scattering experiments. This corresponds to the following approximations in the 
integrand expression in the formula H51|l 

where = hk/m is the group velocity of the center of the wave packet coinciding with the velocity of classical particles. 

Let us remind briefly results of the time-dependent collision theory in the case of the short-range potential when 
the asymptotic form ((SJ of the stationary wave function can be used for analysis of the wave packet evolution 

where 9q is the angle between the vectors q and qr/r. 

Now one can use the expansions H52|l and to find the following result for the function ^'g(r, t) 

vl/,,(f, t) ~ {e'^^^^G(|f - fo - vt\) + f{0,) — Gi\4 - ro ~ vt\)}e-'^''«e~''''''''/^"\ (54) 

V rv 

Fig. 5 shows the sketch of distribution of the probability density corresponding to the wave packet (|54|l in some 
moment t. It demonstrates two essential results which represents actually the basis for use the quantum mechanical 
stationary scattering theory for description of the collisions between real particles • Firstly, the overlapping of the 
fluxes corresponding to the incident ( first term in the formula (|54|l l and scattering particles is essential only in the 
above-mentioned interference zone with the angular width 9int = a/r and they can be considered separately out of 
this domain. Besides, the scattering flux is localized in the spherical layer with the average radius r ~ |ro -f vt\ and 
width ~ a. The angular distribution of the scattering particle in the limits of this layer is completely defined by the 
scattering am; 




FIG. 5: Sketch of distribution of the probability density corresponding to the wave packet l|54^ in some moment t 

The expansions (|52|l can be used in the integral l|53|l in the case of the integrand without singularities in the range 
of the variable variation . This condition doesn't satisfied for the asymptotic form Q in the case of the Coulomb field 
because the Rutherford amplitude includes unintegrable singularity. Let us show, however, that the representation of 
the wave packet analogous to the formula (|54|l is justified also for the Coulomb problem if the expansion H54|l is built 
on the basis nonasymptotic representation (jl^ll for the confluent hypergeometric function: 
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Zq = qr — qr. (55) 

The functions 6*1^2 are rather smooth and integrable. One can use the expansion H52|) for their arguments if the 
foUowing condition is satisfied in the region < 1 of the most essential variation of these functions 

k 

Zk > p{r- - r) ~ prOp. (56) 

If the spread of the wave packet is neglected, the value \p\ can be estimated as |p| ~ kOp ~ ka/r ( is the characteristic 
linear size of the wave packet localization in space) and the condition H56|l leads to the inequality 

02 ^ 1 > («)2^ ^ < 1. (57) 
kr r r 

It coincides with the above mentioned estimation (|49|l considered on the basis of the qualitative analysis. 

In the result the functions 01.2 in the formula H55|) can be removed out of the integral with the arguments corre- 
sponding to the center of the wave packet and it leads to the expression 

^~^{r,t) = [[zkr^>'-^^GMk,^Zk)G{\r^ - fo - vt\) + 

{zk^^'^!^M2{^^k.^Zk)G{\T- To - vt\)]e-''''°e~^'^''"'/'"^; 

Zk = kr — kr. (58) 

It means that the scattering process in the Coulomb field can be considered on the basis of the stationary theory as 
it takes place in the case of the short-range potential. Besides, the incident and scattered wave packets are extending 
in the space separately excluding unessential domain of their overlapping. 

VI. CALCULATION OF THE CHARGE CARRIER MOBILITY IN THE EXTRINSIC 

SEMICONDUCTORS 

It is important to consider the concrete physical system where the described peculiarities of the scattering process 
in the Coulomb field can be appeared for some observed characteristics. Accordingly to the estimation H49|) . it is 
possible if the following inequality is fulfilled 

^ > 1. (59) 
r 

Here is the screening radius of the Coulomb potential in a medium and it depends on the screening mechanism 
in the system. The value r is defined by the distance between the scattering center and detector or by the average 
distance between two subsequent collisions if the scattering operator H30|l is used for the description of kinetic processes 
in the system. 

In the present paper the nonasymptotic scattering theory will be used for analysis of the charge carrier mobility 
in the extrinsic semiconductors for low temperature. In this case concentration of the impurity centers defines both 
the type of the carriers and their concentration and also the main contribution to the resistance of the semiconductor 
|l7j| . The problem was recently analyzed in detail in the paper Q and results of the various phenomenological models 
for regularization of the Rutherford cross-section were compared with the experimental data [7|. It was shown that 
the wide used models of Brooks- Herring 5], and Conwell-Weisskopf Jil don't describe completely the experimental 
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dependence of the mobility on the temperature and impurity concentration. The authors of the paper fitted the 
experimental data essentially better by means of an additional phenomenological parameter with the physical meaning 
of the characteristic time of the collision. It seems to us that such parameter takes into account partly the influence 
of the "near" zone (see Sec. 2) on the formation of the scattered flux. So, the regularization of the scattering problem 
in the Coulomb field is of interest not only as the methodical problem but also as the applied one. 

Let us consider the extrinsic semiconductor with the concentrations of the donors ni and acceptors 712 in the charge 
states Zie and correspondingly (in the most of real structures the impurities with the charge \Zi,2 \ = 1 are mainly 
important ), e is the absolute value of the electron charge. 

In general case the value is defined by both the thermally excited carriers and the carriers conditioned by the 
impurities. The semiconductors with the wide forbidden zone were analyzed in the paper and the value rig can be 
estimated as 

Ue — ZiTLi — Z2n2 = n, 

for the considered low temperature. 

Let us introduce also another parameter which is more spread in the semiconductor physics: K is the compensation 
and is usually a quite small value if ~ 0.1 Q 

n2 ^ — TTT^, K= — . (60) 



Zi — KZ2 Z\ — KZ2 rii 

It is well known that the Coulomb potential screening in semiconductors is defined by several factors. From 
one side, there is the static dielectric constant e conditioned by the electrons from the valency zone which doesn't 
change the long-range character of the potential. From the other side, the Debye screening of the potential by free 
electrons (or holes) leads to its cut off on the distance ^3 



where is the Boltzmann constant; T is the crystal temperature; n^, is the concentration of free charge carriers 
(electrons in the conductivity zone for n-type semiconductors or holes in the valency zone for p-type semiconductors). 

The average distance r between scattering centers and the characteristic wave vector for the carriers in the formula 
(|55)l can be estimated as 



, V2m*E y/3m*kBT 
k^—^^ , 

with m* as the carrier effective mass. 

In the result the condition (|59() leads to the following inequality 

(3em-)V^(fc^T)3/^ > 1, (62) 

which is fulfilled in the entire range of the density and temperature considered in . 

In the most applications the theoretical estimation of the carrier mobility is based on the approximation of relaxation 
time r and the Maxwell velocity distribution. It leads to the following formula (n-type semiconductors are considered 
for the definiteness) 

/X = ^ < r > (63) 

TO* 

< r [ / ij3/2g-B/fc«T]-i / ^~ r{E)E^/^e-^/^^^. (64) 



Here the relaxation time is supposed to be averaged on the energy of carriers with the Maxwell distribution. 
It is known To] that if the several mechanisms of scattering take place (in our case there are scattering by donors 
and acceptors), the more accurate result the additional averaging on the types of scattering centers should be fulfilled: 
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TiAE) = ^ , (65) 

where the indexes 1,2 correspond to the scattering by donors and acceptors; atri,2 is the transport cross-section for 
the cases of the attraction and repulsion. In accordance with Sec. 3 these values are defined by the formulas : 



(ytrl.2 = -^^2 1^(1 ± *^1^2, 1, iz)Yzdz. (66) 

We use here the more accurate formula than Ea. H31|l because in this case the condition /cr ^ 1 can not be fulfilled. 
The parameters of interaction between carriers and scattering centers in the considered cases are the following 

t.l,2{E) - — 

env 

and the static dielectric constant of the crystal is taken into account. 

Accordingly to the formulae l|66(l the transport cross section depends on the potential charge as distinct of its calcu- 
lation with the Rutherford cross section. The similar effect ("phase shift") is well known for extrinsic semiconductors 
and considers usually by means of the Fridel sum rule p^ . Indefinite parameter r is included in Ea. H66|l . If the value 
/I is calculated by the totally microscopic way it should be averaged on the space distribution of the impurities in 
the sample. It is equivalent to the integration of the expression (|63|l by r taking into account Ea. l|66(l . However, the 
transport cross-section has the smooth logarithmic behavior on r which can substituted in Ea. H66|l as the average 

distance between the impurities with the considered accuracy. Then the value r = O.Sn^ can be used in Ea. H66f) 
analogously to the both models 5] , and fil . 

It is convenient to define the auxiliary value so, that 



^2 



Then nonasymptotic calculation leads to 



r-2fcri,2 

eT"«i-^ / |C/(l±ia,2,l,iz)|2zdz, (67) 



with the values 



^ (Zi - KZ2Y'^ ^ {Z^-KZ2Y'^ 

''^ 2ni/3 ' 2(nif)i/3 ' 

which are defined by the half of the average distance between the donors and acceptors correspondingly. 
In the result the following expression for the carrier mobility can be obtained: 

^ 3^3/2e3^*i/2„ 7^ Zlo[^,{x)^KZlo[^^{xf''- ^^^^ 

The integrals over energies can be estimated by the standard way 6] : the smoothly changing functions can be taken 
out of the integrals with argument x = 3 when the energy distribution function has the maximum value. It leads to 
the following analytical expression for the mobility: 



2'^l''e^{Z^^KZ2){kBT) 



3/2 



^3/2e3m*i/2„(z2CT^,i(3A:BT) + KZlu'^^^^iikBT)) 



(69) 
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In the case Zi = Z2 = 1 it transforms as follows 

27/V(i-x)(fcBr)3/2 



We can compare it with the analogous formula in the framework of the Conwell-Weisskopf model 



(70) 



The results of calculation by means of Eas. H70|l and their comparison with the Conwell-Weisskopf model results are 
shown in Fig. 6. The same figure shows that the dependence of the mobility on the temperature and compensation K 
in our consideration differ essentially on the results of the Conwell-Weisskopf model based on the Rutherford cross 
section with the phenomenological regularization. In principle, such distinction can be discovered in some experiments. 




FIG. 6: Comparison of the mobilities calculated with the nonasymptotic transport cross-section (solid line) and in the frame- 
work of the Conwell-Weisskopf model (dashed line) (6a - dependence on the impurity concentration; 6b - on the compensation; 
6c - on the temperature). The following parameters were used T = 78, e = 10,m — 0.2mo,K — 0.15. 
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